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^ Abstract 

We suggest a simple modification of the maximal mixing scenario (with S3 
permutation symmetry) for three light neutrinos. Our neutrino mass ma- 
trix has smaller permutation symmetry S2 iyn <-> v e ), and is consistent with 
all neutrino experiments except the 37 C1 experiment. The resulting mass 
eigenvalues for three neutrinos are mi ~ (2.55 — 1.27) x 10~ 3 eV, 7712,3 ~ 
(0.71 - 1.43) eV for Am 2 LSND = 0.5 - 2.0 eV 2 . Then these light neutri- 
nos can account for ~ (2.4 — 4.8)% (6.2 — 12.4%) of the dark matter for 
h = 0.8 (0.5). Our model predicts the Un — * v T oscillation probability in the 
range sensitive to the future experiments such as CHORUS and NOMAD. 
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The minimal standard model (MSM) has been highly successful in describing interactions 
among elementary particles from low energy up to ~ 100 GeV. The only possible exception 
may be various types of neutrino oscillation experiments. There have been positive indica- 
tions from large scale experiments for solar and atmospheric neutrinos that a certain amount 
of mixing between neutrino species may be present jl]]. The recent report from the LSND 
experiment at the laboratory scale provides us with another hint of such a possible neutrino 
mixing [0] . Since neutrinos in the MSM are exactly massless, there can be no mixing among 
them, and it is impossible to accommodate such neutrino mixing data in the framework of 
the MSM. This situation is rather encouraging, since it is at present the only place where 
we can grasp a hint of new physics beyond the MSM. 

In view of this, it is quite interesting to speculate what type of neutrino mass matrix can 
fit all the data from the various types of neutrino oscillation experiments. It is our purpose 
to present one such mass matrix in this work. Most analyses on the neutrino oscillation 
assume two neutrinos oscillating with one mass difference parameter, Am 2 . However, the 
LSND experiment and the atmospheric and solar neutrino data hint at least two mass 
difference parameters, requiring oscillations among at least three neutrinos. For oscillations 
with three neutrinos, we have two mass differences, three real angles and one phase. In 
order to simplify the analysis, a certain ansatz for the mass matrix is required. In this 
vein, we first briefly discuss the maximal mixing scenario for the neutrino sector. We then 
present our ansatz for the neutrino mass matrix as well as the numerical analyses to fit the 
atmospheric, LSND, and solar neutrino data from GALLEX and SAGE. In this work, we 
consider oscillations among three neutrinos only, v a — > vp. 

One of the popular ansatz for the neutrino mass matrix is the maximal mixing one 
(equivalent to a cyclic permutation symmetry among three generations) || : 




M maximal = | b* a b I , (1) 
with the mixing matrix U given by 



UJl UJl UJl 
UJi UJ 2 UJ 3 
UJi UJ 3 UJ 2 



U maximal ~ \ ^1 &2 ^>3 | • (2) 



Here, ^2,3 are three complex roots of uj 3 = 1 with uj\ — 1. This ansatz was originally 
proposed in the neutrino sector || , and extended to the quark sector with a partial success 
in explaining the quark masses and the Kobayashi-Maskawa matrix elements [|J. 

The maximal mixing scenario has many interesting features || . For example, the survival 
probability for a neutrino is independent of its flavor. The v e survival probability has 
two plateaus, 5/9 in the intermediate step, and 1/3 for -| 3> (A \ , through vacuum 

oscillations. With vacuum oscillation, one cannot explain solar neutrino data from 37 C1 
and the Ga data simultaneously. Thus one must make a choice between the 37 CI and Ga 
experiments. Here, we choose to interpret the Ga data (GALLEX and SAGE experiments) 
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through vacuum oscillation and disregard the 37 CI data [], following Ref. 0. The Ga data 
requires v e survival probability of ~ 5/9, which implies ~ for the solar neutrino is smaller 
than jjr — k — — . The atmospheric neutrino data and km range laboratory experiments require 



%3' 

,2 



another Am|. Hence, two mass difference scales have been all used, and there is none left 
for a new scale suggested by the LSND data around Am 2 LSND ~ 0(1) eV 2 with a mixing 
angle ~ (a few) x 10~ 3 . The only possibility to explain both the mass shifts at LSND 
point and at atmospheric data points in the maximal mixing scenario is that there are two 
thresholds corresponding to a larger Am 2 at ~ 0(1) eV 2 and a smaller Am 2 at around 10 -2 
eV 2 . In this case, the v e survival probability for the solar neutrino problem is 1/3, and is 
too small to accommodate the Ga data. Therefore, although the qualitative features of the 
maximal mixing scenario is encouraging, it is not viable if the LSND data is confirmed in 
the future. Another way to see this is as follows : the maximal mixing scenario predicts 
the transition probability for — * u e to be 4/9 in the range of the LSND experiment, 
which clearly contradicts the reported transition probability, (a few) xl(R 3 . Furthermore, 
the best x 2 fit to the atmospheric and the solar neutrino data indicates that the masses of 
three light neutrinos are m 3 ~ (85 ± 10) meV, and m 12 < 3/ieV ||, which are too light to 
be cosmologically interesting as a hot dark matter component of the missing mass of the 
universe. 

Therefore, we make an ansatz for the neutrino mass matrix which is a simple modification 
of the maximal mixing one, (1), and study its consequences in this work. We assume that 
neutrinos are Dirac particles so that the lepton number is to be conserved in our model. 
Then, each left-handed neutrino (v l L ) is accompanied by the right-handed partner (v l R ) which 
is sterile under electroweak interactions. 

Note that any 3x3 matrix Mij can be decomposed as M = X — iY with both X and 
Y hermitian. Also any hermitian matrix X can be written as X = S + iA, where S (A) 
is a real (anti) symmetric matrix. Finally, the symmetric matrix S can be decomposed as 
the trace part proportional to 5ij and the traceless symmetric matrix. One can combine the 
trace part of the symmetric mass matrix 5* and the real antisymmetric part A in order to 
get the neutrino mass matrix, 



M = fj, 



/ 1 ic id 

—ic 1 ib ] , (3) 
y —id —ib 1 



where /i is the mass scale, b, c and d are all real. We have chosen a basis in which the charged 
lepton mass matrix is diagonal. Note that the diagonal terms are still universal, and only the 
off-diagonal elements are modified from the maximal mixing one, (1). Note that this ansatz 
becomes the maximal mixing one if b = c = —d (i.e., if there is a permutation symmetry 
among three generations). This form for the mass matrix is sufficiently simple but rich 



The matter oscillation effect, i.e., the MSW mechanism, has also been suggested to interprete 
this Homestake experiment in Ref. ||. In this case, the relevant Am 2 is around ~ 10~ 4 eV 2 . Since 
we have only two Am 2 around ~ 0(1) eV 2 and ~ 10~ 2 eV 2 , the MSW mechanism is irrelvant to 
our study in this work. 
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enough to give nontrivial analytic formulae for the survival and transition probabilities for 
three neutrinos. If there is a solution when two of the off-diagonal elements are the same 
(say, b = d for example), then the permutation symmetry among three generations (S 3 ) 
in the original maximal mixing case (1) breaks down to S 2 . This would imply that our 
mass matrix ansatz depends on three real parameters and, thus it is one of the simplest 
modifications to the maximal mixing ansatz, which can accommodate LSND, atmospheric 
and solar neutrino data. In fact, this is the case (with b 2 = d 2 ) as discussed in the following. 
It breaks the original S3 possessed by (1) into S 2 in a particular way. There may be several 
other ways to break this permutation symmetry, which will be considered elsewhere. There 
may be some underlying dynamical reasons for the above form of the mass matrix, but we 
take it as a simple phenomenological ansatz for the moment. 
Three eigenvalues of the mass matrix (3) are 



mi = fx, 



m 2 , 3 = /i(l± AO, 



(4) 



where N = (b 2 + c 2 + d 2 ) 1 ^ 2 . The corresponding eigenvectors form the mixing matrix U 
which relates the weak eigenstate v a to the mass eigenstates as v a = Y^iUaiVi. The indices 
a = e,u,r label the flavor eigenstates, and i = 1, 2, 3 label the mass eigenstates of three 
neutrinos. 

Then one can easily verify that 



P(v„ 



P(v„ 



1 I 

N 4 



V e ) = 7 

; N 4 



(A 21 + A 3 i) b 2 (N 2 ~b 2 ) + - A 32 (N 2 - b 2 ) 2 
1 



P(V, -!/„) = 1 




Vr) = w 



(A 21 + A 31 )b 2 d 2 + - A 32 (c 2 N 2 - b 2 d 2 ) 



(A 21 + A 31 )d\N 2 - d 2 ) 

2 J + 2{c 2 + d 2 ) 2 

(A 21 + A 31 )cV + A 32 - 



(5) 
(6) 
(7) 



2(c 2 + d 2 y 



(8) 



where 



9 ( 1.27 L Am 2 / 
A,, =2 sin 2 ^ 



(9) 



with Am 2 - = m 2 
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0, there exist only two 



in eV 2 and L/E in km/GeV. Since E Am?- 
independent mass difference parameters. 

Note that the heights of the plateaus for the v e survival probability are functions of 
b 2 /N 2 only, even in the presence of nonvanishing c and d. Since Piy^ — > u e ) depends on 
two mass differences, one can identify the mass difference Am 2 LSND ~ 0(1) eV 2 either as 
Amjj or as Am 2 2 . The other mass difference is taken to be 0.72 x 10~ 2 eV 2 in order to solve 
the atmospheric neutrino problem. Thus, there is a reasonable hierarchy between two mass 
differences. In the following, we discuss two possibilities separately. 
(I) : Am 2 x = Am 2 LSND ~ 0(1) eV 2 and Am 2 32 = 0.72 x 10~ 2 eV 2 : 
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In this case, the heights of the intermediate plateaus for the v e and survival probabil- 
ities are given by 



b 2 ( b 2 \ 

P{v e ->v e ) = l-2 — (l-^J, (10) 

d 2 ( d 2 \ 

Piy, - uj = 1 - 2^ ^1 - _ j , (11) 

and the transition probabilities can be approximated as 

b 2 d 2 2 /l.27LAm 2 \ , 
P(i/„ - v.) = 4— sin 2 ^ , (12) 

v c 2 d 2 . 2 /l.27LAm 2 \ 
P(u, - „ T ) = 4— sin 2 f . (13) 

The transition probability is often described in terms of two parameters, Am 2 and 9 a p for 
which 

P[v a - ^) = sin 2 2^ sin 2 ^- 27 ^^ . (W) 

Therefore, we can make the following identifications : 

Ab 2 d 2 

sin 2 2^ = ^, (15) 
4r 2 d 2 

sin22 V = ^r, (16) 

in our model (for Am^ >> Am§ 2 ). 

Experimental results for the neutrino oscillations are shown in the (sin 2 29, Am 2 ) plane. 
In the plot presented by the LSND group, there are small regions in this plane which indicates 
a possible transition of — > 9 e . This region is not ruled out by other laboratory searches 
such as BNL E776 0, KARMEN §, BUGEY § and others 0- 0. For each possible 



Am 2 LSND , we show the possible value(s) of sin 2 29 e ^ in Table 1. For the same Am| 



LSNDi 



there is an upper bound on sin 2 MT from FNAL E531 @, CHARM II |3fl and CDHSW 
TM, and we also list these numbers in Table 1. 



For each Am^ = /S.m 2 LSND given in Table 1, one can solve Eq. (4) to get the neutrino 
masses. For example, Am| x = 6 eV 2 leads to 

m 1 = 7.35 x 1(T 4 eV, 

m 2 « -m 3 « 2.45 eV, (17) 

with Si|m Wi | = 4.9 eV. (The negative can be remedied by a chiral transformation of 
1/3 field.) For other values of Am 2 1; we show the resulting neutrino masses from our mass 
matrix ansatz (3) in the fourth column of Table 1. 

These light neutrinos can contribute to the missing mass of the universe (the hot dark 



matter) in amount of 15 
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Clh 2 = 7.83 x ID- -ML- fe) , (18) 
g*s{T D ) \eVJ 

where g* s (T D ) = 10.75 and g e ff = (3g)/4 = 3/2 are the effective degrees of freedom 
contributing to the entropy density s and to the ratio Y = n/s, n being the num- 
ber density, respectively. The parameter h is related to the Hubble constant Ho as 
H = lOO/i km sec -1 Mpc -1 [HJ. So, for the solution (17), three light neutrinos can consti- 



tute 8.3 % (21.4 %) of the missing mass of the universe for h = 0.8 (0.5), which is again 
cosmologically interesting. (Here, we have assumed that three sterile right-handed neutri- 
nos decouple much earlier than the left-handed neutrinos, and that they don't affect the 
results of the standard cosmology.) The results for other values of Am^ are listed in the 
last column of Table 1. 

When we determine b 2 and d 2 , it is important to satisfy all the constraints shown in Ta- 
ble 1. One might try to perform the \ 2 fit to the available data on the neutrino oscillations. 
Instead, we choose to scan d 2 /N 2 for each Am 2 in the first column of Table 1. For each 
d 2 /N 2 , the parameter b 2 is determined by the mixing angle given by the LSND experiments, 
and c 2 /N 2 = (1 — b 2 /N 2 — d 2 /N 2 ). Then, we require that the resulting sin 2 2# MT satisfy the 
upper limit given in the third column of Table 1. We also calculate the survival probabilities 
for u e and at the intermediate level (the laboratory and the km range scale) and require 
them to be larger than 0.95 in order to satisfy the null results in various types of disappear- 
ance experiments for the and v e beams. For Am 2 = 6 eV 2 or larger (the first and the 
second rows), the resulting d 2 /N 2 0.995 — 1.00, which corresponds to almost no disappear- 
ance of i/p for all ranges of L and E. Thus, we reject d 2 /N 2 around 1. For Am^ < 2 eV 2 , 
the allowed ranges for d 2 /N 2 are typically around 0.010-0.020. The corresponding b 2 /N 2, s 
are also in the same range as d 2 /N 2 . Thus, as discussed in the following, we can accommo- 
date the laboratory scale and the large scale neutrino experiments, by choosing small (but 
nonvanishing) b 2 and d 2 . 

In particular, there is a small region in which b 2 = d 2 gives acceptable fits to all available 
data on neutrino oscillation experiments except for the 37 C1 data. This is quite interesting, 
since it corresponds to residual permutation symmetry (S 2 ) between e and \i in (3) with three 
real parameters. In other words, we have \M efl \ = |M MT |. Our matrix with b 2 = d 2 breaks 
the original symmetry of the maximal mixing one (S3) into S2, and thus may be regarded 
as one of the simplest modifications to the maximal mixing one, (1). In the following, we 
demonstrate that the ansatz (3) can reasonably fit all the data with b 2 /N 2 = d 2 /N 2 = 0.015 
for Affijj = 2 eV 2 , except for the HOMESTAKE data, with a reasonable accuracy at present. 

In Figure 1, we show the resulting survival probability for u e in the solid curve Q along 
with various types of neutrino oscillation data, the v e disappearance experiments at reactors 
||- [[II] and the solar neutrino experiments ||16|| - ||19|| . Using three neutrino mixing with four 



parameters, we get two step survival probability for v e — > v e . The plateau for the large L is 
about 0.49, a little bit lower than 5/9 used before in the maximal mixing case. Thus, the 
solar neutrino deficit is solved in terms of vacuum oscillations. The intermediate plateau of 



2 In order to average the oscillation probabilities, we adopt the prescription by Harrison et al. 
which amounts to replacing cos(:c/2) by smx/x. 
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~ 0.97 for v e is the prediction of our specific form of the mass matrix, and is consistent with 
all the existing data. The KRASNOYARSK data has a relatively large error bar, and our 
curve is within two a of the data. 

In Figure 2, the survival probability for using our ansatz is shown in the solid curve, 



along with the disappearance experiment data ||21|| . The survival probability for z/ M 
at large L is 0.48, and the intermediate plateau has a height of 0.97. The curve agree with 
the available data quite well. 

For the atmospheric neutrino data, we show the data point for the so-called R defined 

by 



along with our prediction 



Data /i 

= (N,/N e ) MC ' Uyj 



R = P ee + rP, e (20) 



in Table 2, where r is the incident (/Ve) ratio. From Table 2, we observe that most of the pre- 
dicted R values are consistent with all the atmospheric neutrino data from KAMIOKANDE, 



1MB, and others |p2]| - pq| , considering some data points have large errors. 

For these numbers with Am^ = 2 eV 2 , we predict sin 2 2$^ m 6 x 10~ 2 , which is 
just below the current upper limit, 8 x 10~ 2 . This range of sin 2 26 may be probed at 
CHORUS, NOMAD, FNAL P803, CERN/ICARUS and FNAL/SOUDAN2 |7fl. It would 
be interesting to test our predictions for the — > v T oscillation in the future. Similar results 
can be drawn for other values of Am^. Thus, our mass matrix ansatz (3) not only describes 
all the available data on neutrino oscillations, but also predicts the mixing angle for — > v T 
in an interesting range which lies within sensitivity of the near-future experiments. 

Let us briefly discuss the second case : (II) Am 2 x = 0.72 x 10~ 2 eV 2 , and Am 2 2 = 
Am 2 LSND . In this case, it is easy to verify that there is no solution for b 2 /N 2 and d 2 /N 2 
which satisfy the constraints from the laboratory scale experiments from BUGEY, BNL 
E766, and those in the second and the third columns of Table 1. So, our mass matrix ansatz 
prefer the solution (I) for which mi is smallest around 10~ 3 — 10~ 4 eV, and the other two 
are nearly degenerate with 1712 ~ w O(l) eV. 

In summary, the neutrino mass matrix ansatz (3) with four real parameters is one of the 
simple modifications to the maximal mixing ansatz (that fails to fit the new LSND data) 
which can fit various types of neutrino experiments except for the 37 CI solar neutrino data 
17 1 . In particular, there are solutions with b 2 = d 2 with residual permutation symmetry 



among two generations {y^ «-> u e ). In this sense, our mass matrix ansatz could be regarded 
as one of the simplest modifications to the maximal mixing ansatz. The resulting light 
neutrinos have masses (17), and thus they can constitute about 2.4-4.8 % (6.2-12.4 %) of 
the missing mass of the universe for h = 0.8(0.5), and thus cosmologically interesting unlike 
the maximal mixing scenario. In our model, the transition probability for — > v T is close 
to the current upper limit, depending on Am 2 LSND as shown in the third column of Table 1. 
Since it lies within the reach of various future experiments such as CHORUS, NOMAD, 
etc., the observation of — v T oscillation would constitute a definite test of our mass matrix 
ansatz along with the confirmation of the LSND data. 
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Note Added in Proof 



In Eqs. (5)-(8), we did not show terms involving = sin (^——^—^j, since these terms 
vanish under taking averages, or they are irrelvant to our study. See Ref. [5] for more details. 
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FIGURES 



FIG. 1. The survival probabilities P(y e — > v e ) using our ansatz (3) for the case (I), along with 
the reactor experiment data from KARMEN, ILL/GOSGEN, BUGEY, KRASNOYARSK, and the 
solar neutrino data from KAMIOKA, HOMESTAKE, SAGE and GALLEX. 

FIG. 2. The survival probabilities of using our ansatz (3) for the case (I), along with the 
accelerator experiments from CDHS-SPS and CHARM-PS. 
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TABLES 



TABLE I. The allowed regions for Am 2 and sin 2 29 e ^ consistent with the LSND as well as 
BNL E766, and the corresponding upper limit for sin 2 29 ^ T from FNAL E531 and CDHSW. The 
fourth column is the predicted neutrino masses by our mass matrix ansatz (3). The last column 
shows contributions of three light neutrinos to the missing mass of the universe for h = 0.8 (0.5). 
See the text for details. 



Am 2 (eV 2 ) 


sin 2 29 efl 


sin 2 29^ 


(mi, m 2, m z) m e V 


n (%) 


20 


~ 3 x 10" 3 


< 4 x 10" 3 


(4.02 x 10" 4 , 4.47, -4.47) 


15.1% (38.9%) 


6 


2 x 10~ 3 


< 2 x 10~ 2 


(7.35 x 10~ 4 , 2.45, -2.45) 


8.3% (21.4%) 


2 


(1 ~ 2) X 10~ 3 


< 8 x 10~ 2 


(1.27 x 10 -3 , 1.43, -1.43) 


4.8% (12.4%) 


1 


(2 ~ 6) x 10~ 3 


< 0.1 


(1.80 x 10 -3 , 1.0, -1.0) 


3.4% (8.8%) 


0.5 


(0.7 ~ 2) x 10~ 3 


< 0.2 


(2.55 x 10 -3 , 0.71, -0.71) 


2.4% (6.2%) 



TABLE II. The atmospheric neutrino data R for various L/E along with our predictions for 
Am 2 ,! = 2 eV 2 , Am| 2 = 0.72 x 10~ 2 eV 2 and 6 2 /iV 2 = d 2 /N 2 = 0.015. We show the r = (fi/e) incident 
values for each data point also. 



Experiments 


r 


L/E (km/GeV) 


Measured 


Prediction 


KAMIOKA (2| 


4.5/1 


5 


i 97 +0.61 
■ L ' z '-0.38 


0.99 


(Multi-GeV) 


3.2/1 


10 


6S+ - 21 
U - D,5 -0.16 


0.97 




2.2/1 


100 


o-5i±g:ii 


0.41 




3.2/1 


1000 


0.46l°i« 


0.31 




4.5/1 


2000 


0.28i°; ? 


0.22 


KAMIOKA || 


2.1/1 


80 


0.59 ±0.10 


0.50 


(Sub-GeV) 


2.1/1 


12800 


0.62 ±0.10 


0.48 


1MB H 


2.1/1 


1000 


0.54 ±0.13 


0.47 


FREJUS || 


2.1/1 


500 


0.87 ±0.18 


0.47 


NUSEX H 


2.1/1 


500 


0.99 ±0.32 


0.47 


SOUDAN H 


2.1/1 


1000 


0.69 ±0.21 


0.47 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503436v4 
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